Incomplete equilibrium in long-range interacting systems 
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We use a Hamiltonian dynamics to discuss the statistical mechanics of long-lasting quasi- 
stationary states particularly relevant for long-range interacting systems. Despite the presence of an 
anomalous single-particle velocity distribution, we find that the Central Limit Theorem implies the 
Boltzmann expression in Gibbs' F-space. We identify the nonequilibrium sub-manifold of F-space 
characterizing the anomalous behavior and show that by restricting the Boltzmann- Gibbs approach 
to this sub-manifold we obtain the statistical mechanics of the quasi-stationary states. 
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In comparison with its equilibrium counterpart, 
nonequilibrium statistical mechanics docs not rely on uni- 
versal notions, like the ensembles ones, through which 
one can handle large classes of physical systems (IJ] . In- 
complete (or partial) equilibrium states [2,13 are in this 
respect a remarkable exception, since in these cases con- 
cepts of equilibrium statistical mechanics can be used 
to describe nonequilibrium situations. Incomplete equi- 
librium states arise when different parts of the system 
themselves reach a state of equilibrium long before they 
equilibrate with each other . The classical understand- 
ing on how a system approaches equilibrium is based on 
the short time-scale collisions mechanism which links any 
initial condition to the statistical equilibrium. For long- 
range interacting systems, this picture is not valid any- 
more since the time-scale for microscopic collisions di- 
verges with the range of the interactions. This implies 
that the Boltzmann equation must be substituted with 
other approximations such as the Vlasov or the Balescu- 
Lenard equations 3 , where the interparticle correlations 
are negligible or almost negligible and a nonequilibrium 
initial configuration could stay frozen or almost frozen 
for a very long time. This applies, e.g., to gravitational 
systems, Bose-Einstein condensates and plasma physics 
[3 . Due to the physical relevance of long-range interact- 
ing systems and to the privileged position of incomplete 
equilibrium states in nonequilibrium statistical mechan- 
ics, it is important to investigate whether the notion of 
incomplete equilibrium plays an important role in under- 
standing the nonequilibrium properties of these systems. 

Recently we showed @ that nonequilibrium states 
in which the value of macroscopic quantities remains 
stationary or quasi-stationary for reasonably long time 
(quasi-stationary states - QSSs) are important, e.g., for 
experiments, since they appear even when the long-range 
system exchanges energy with a thermal bath (TB). Us- 
ing the same paradigmatic long-range interacting system 
of Ref. 0, the Hamiltonian Mean Field (HMF) model 
0, here we discuss the Gibbs' F-space statistical mechan- 



ics description of the QSSs in a canonical ensemble per- 
spective. We identify the nonequilibrium sub-manifold 
of F-space within which the quasi-stationary dynamics is 
confined and we show that the Boltzmann-Gibbs (BG) 
approach, restricted to this sub-manifold, gives the cor- 
rect statistics of the QSSs. In this respect, the QSSs can 
be interpreted as incomplete equilibrium states j^]. Our 
theoretical framework allows one to calculate, on the ba- 
sis of the empirical detection of the temperature and of 
the value of an order parameter, any other thermody- 
namic quantity such as the energy or the specific heat of 
the system. The possibility of predicting physical quanti- 
ties which characterize the QSSs could be useful, i.e., for 
understanding nonequilibrium features of gravitational 
or plasma structures and it is then of particular interest 
for experimentalists or theorists of long-range interact- 
ing systems. Since the system considered is naturally 
endowed with a microscopic Hamiltonian dynamics, we 
validate step by step our theoretical derivation with a 
priori results obtained from dynamical simulations. Our 
findings also furnish novel significant arg uments to an in- 
tense debate in the literature H, 0, llOl ITl| , that so far 
has been restricted to the single-particle /i-space and to 
the microcanonical ensemble. 

The HMF model can be introduced as a set of M glob- 
ally coupled XK-spins with Hamiltonian 
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where 6i g [0, 27r) are the spin angles and S R their 
angular momenta (velocities). The specific magnetiza- 
tion of the system is tohmf = | X]f=i("^os^i'™^i)l/-^ 
and the temperature T is identified with twice the spe- 
cific kinetic energy. We have thus chmf = 7hmf/2 -I- (1 — 
"^hmf)/2i where chmf = -Ehmf/A/ is the specific energy. 
Direct connections with the problem of disk galaxies [ij 
and free electron lasers experiments 13] have been estab- 
lished for this Hamiltonian. Eq. has also been shown 
to be representative of the class of Hamiltonians on a one- 
dimensional lattice in which the potential is proportional 
to X]f j=i [1 ~ cos{6i — 6j)] /rf^, where 7^ is the lattice 
separation between spins and a < 1 [lj|. Hence, the 



Hamiltonian in Eq. can be considered as an interest- 
ing "paradigm" for long^ange interacting systems (l^. 
The TB introduced in Q is characterized by iV 3> M 
equivalent spins first-neighbors coupled along a chain 
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with 9n+i — Qm+i, and the interaction between HMF 
and TB is given by 
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where e is a coupling constant that modulates the inter- 
action strength between HMF and TB. Each HMF-spin 
is thus in contact with S TB-spins specified as initial 
condition {rs{i) are independent integer random num- 
bers in the interval [M -|- 1, A^]). A "surface- like effect" 
S ~ M'*~^ (0 < 7 < 1) guarantees a consistent thermo- 
dynamic limit . For 6 = HMF and TB are decoupled 
and the setup reproduces the microcanonical dynamics. 
For e 7^ the whole system is at constant energy whereas 
the energy of the HMF model fluctuates. Our numerics 
are obtained with M = 10^, N = M^, S = lO^M-^^^, 
0.005 < e < 0.1 (we use dimensionless units), through 
a velocity- Verlet algorithm assuring a total energy con- 
servation within an error AE/E < 10~^ ^6j. The width 
To of the Maxwellian probability density function (PDF) 
for the initial TB- velocities is a control parameter for the 
bath temperature. For e > we showed that the HMF 
temperature finally converges to the BG equilibrium at 
temperature Tq. 

By setting far-from-equilibrium initial conditions for 
the HMF model, the relaxation to equilibrium typi- 
cally displays stationary or quasi-stationary stages dur- 
ing which the phase functions mHMF, Thmf (and thus 
also chmf) fluctuate around constant or almost constant 
average values . This behavior is particularly interest- 
ing when the life-time of the QSS diverges in the ther- 
modynamic limit 0, H, 0, Hol llTI . This happens if for 
example at t = we set a delta distribution for the an- 
gles (phmf(^) = <5(0) =^ iTT'HMF ~ 1)' ^ uniform dis- 
tribution for the velocities, phmf(0 = 1/2^, I & [—ij], 
with r~ 2.03 (chmf — 0.69) @, and a TB temperature 
To = 0.38. In Fig. [TJi we show that during the QSS, for 
e > 0, the single particle velocity PDF is non-Maxwellian 
and similar to the distribution found in the microcanon- 
ical case (e = 0). 

Given some probability distribution for the initial data, 
a dynamical estimation of phase functions, like e.g. the 
energy E'hmf, can be obtained by recording the phase 
function values at different times in a single orbit and 
averaging over different realizations of the initial condi- 
tions. To understand the connection between the anoma- 
lous PDF in /x-space and the F-space statistics, we start 
by measuring the PDF of the sum of the velocities of L 
particles, Phmf (Fig- Ct*)- Such a distribution tends very 
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FIG. 1: QSS for M = 10=* and To = 0.38. For e = 0.005 we 
observe the average values m^^jp ~ 0.015 and Thmf — 0.397 
(chmf — 0.691). (a): Single particle velocity PDF. The solid 
line is ptb(0- (b): PDF of the sum of the velocities of L 
particles. By multiplying the PDFs for different L's by 1/°'^ 
and dividing the velocities / by L"'^, all data collapse fairly 
well on to the line that corresponds to a Gaussian distribution 
with width T — 0.397. (c): Temperature time-evolutions for 
three different subsets of the system during the QSSs. Here 
the TB temperature has been shifted to To = 0.42. 



quickly to the Gaussian form as L increases. In fact, a 
rescaling of I by L^/^ and a multiphcation of Phmf(0 
by the same factor reveal the Central Limit Theorem 
(CLT) data collapse onto the Maxwellian (Gaussian) dis- 
tribution of temperature T = Thmf = 0.397. The fact 
that the CLT applies to the sum of the velocities is a 
strong indication [i^ that in F-space the probability for 
the energy ii^HMF is characterized by the Boltzmann ex- 
pression w(i?HMF)e~'^'*"''/^ {kB = 1), where lu{Eiimf) 
is a density of states. Although this situation resembles 
equilibrium, there are some important differences. For 
example, the anomalous velocity PDF in /i-space implies 
that the joint probability of all particles is not given by 
a mere product of exponentials. The Boltzmann expres- 
sion arises, because of weak enough particle-particle cor- 
relations for a sufficiently large number of particles. 
Below, we directly verify its occurrence. 

Another key observation is that during the QSS the 
HMF does not thermalize with the TB. In Fig. ^ we 
shifted the TB temperature by 10%, setting it to Tq — 
0.42. While this modifies the final HMF equilibrium tem- 
perature, it does not affect Thmf during the QSS. Even 
the subset of S TB-spins in direct contact with the HMF 
model, {Or,{i)}i<s<s,i<i<M , is at Trfg = Tq and does not 
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FIG. 2: (a); Caloric curve of the HMF model for M = 10^ 
Solid line is the BG equilibrium and dashed line is the curve at 
fixed — 0.015. (b) Theoretical calculation of ln[cj(£'HMF)] 
by using Eq. I^J with the equilibrium caloric curve (solid line) 
and with the curve at — 0.015 (dashed line). 



thermalize with the HMF temperature. The energy fluc- 
tuations are nevertheless significantly larger than those 
due to the algorithm precision (Ai^HMp/^-HMF — 10~^ 
for M — 10"^), distinguishing the canonical QSSs from 
the microcanonical ones. 

We now address the main result of the paper, which 
is central to the discussion of the appropriate statistical 
mechanics approach for quasi-stationary nonequilibrium 
states in long-range systems and to the debate in [slI^ITol 
[Tll |. According to BG, the equilibrium PDF of the energy 
E for a system in contact with a TB at temperature 
T is pbg{E) = uj{E)e-'^/^ /Z, where Z is the partition 
function. Since the Hamiltonian simulations consent an 
empirical estimation of this PDF, it is possible to verify 
Pbg(-E') on dynamical basis ^Jj. From the analytically 
known solution of the HMF model 0, 0| one obtains the 
BG equilibrium caloric curve of the system T[E) (full 
line in Fig. The integration of the thermodynamic 

relation d\nuj{E)/dE = 1/T{E), 



HioiE)] - ln[c^(£;o)] 



dE' 



(4) 



furnishes an analytical evaluation of a;(i?) (full line in Fig. 
|2b) and hence of pbg(-E) [3. In Fig. |3i we show that, as 
expected, pbg(-E'hmf) and the result of the simulations 
at equilibrium, p(£'hmf)i do coincide. A linear regression 
of ln[p(£'HMF)/w(£'HMF)] vs -Ehmf with a coefficient R — 
—0.99997 gives a direct evidence of the Boltzmann factor 
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FIG. 3: For iV = 10^, M = 10^ and e = 0.005 comparison 
between the dynamically recorded p{Eiimf) (empty circles) 
andpBG(SHMF) (dashed lines). 



(Fig. ^p). Moreover, the inverse of the slope coefficient 
agrees with the dynamical T = 2fcHMF within an error 
AT/T = 0.3%. 

With respect to the QSS, it is interesting to ask 
H) ^1 if there exist a statistical mechanics ap- 
proach that, equivalently to the BG equilibrium one, can 
reproduce the dynamically observed p{Eumf)- We first 
notice that the anomalous dynamical behavior during the 
QSS is due to the fact that the system, instead of explor- 
ing the overwhelming majority of P-space microstates, is 
trapped ^1% in regions characterized by almost constant 
nonequilibrium values of the order parameter m. Let 
(to) be the average value around which to fluctuates and 
uj(rn) (E) the sub-manifold of P-space which corresponds 
to this dynamical behavior. The assumption of weak 
correlations among particles, consistent with the previ- 
ous argument based on the CLT and with the Vlasov 
and Balescu-Lenard kinetic pictures , suggests that the 
Lebesgue measure of (E) is non-zero. We then ex- 
pect p{E) = pBG,{m)iE) = o.(™)(^)e-^/^/Z [H. Hav- 
ing assumed this, a saddle point calculation at fixed 
TO = (to) (large deviation formulation of the canonical 
ensemble |l8j at m = (m)) implies that T in the previous 
expression satisfies the fundamental thermodynamic rela- 
tion dliiLL!(^„i^{E)/dE\E={E) = ^/T, where (E) is the av- 
erage value of the energy during the QSS. Hence, uj(„i) {E) 
can be calculated by replacing the equilibrium caloric 
curve T[E) with the caloric curve at constant m = (to), 
T{m){E), and by performing the approach of Eq. (0}. 
The validity of this strategy, and in particular of Eq. 
(4) for the QSS, is further established by the compari- 
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son with the dynamical simulations results. Specifically, 
we show below that T corresponds to twice the specific 
kinetic energy of the HMF. 

The HMF caloric curve at fixed mjjMF = ("^hmf) is 
given, for all (tohmf) G [0, 1], by the straight line 

?HMF,{m>(-E'HMF) = 2^^SMI _ (1 _ (mHMF)^) (5) 

(e.g., dashed line in Fig. |3i for the QSS described in 
Fig. The integration of the inverse of THMF,(m) 

gives W(to)(£'hmf) (dashed fine in Fig. The 
leading behavior of ln[ti;^m> (^hmf)] is proportional to 
Af. This implies that only an exponential probabil- 
ity for the microstates can balance this M-dependency, 
to yield an intensive temperature through the relation 
d\'auj(^rn)(E-mAF)/dE-aM¥- In Fig- lit it is shown that 
p(^'HMf) observed during the QSS at constant {m?) ~ 
0.015 and (2fcHMF} — 0.397 agrees with PBG,(m) (^hmf)- 
Again, a linear regression of ln[p(i?HMF)/w(m)(i?HMF)] 
versus i?HMF with a coefficient R ~ —0.99997 confirms 
the Boltzmann factor for the energy PDF during the QSS 
(Fig. The inverse of the slope coefficient T concurs 
with (2fcHMF) within an error AT/T = 0.5%. We checked 
that a replacement of the limit a — > in the exponential 
Boltzmann factor limQ_>o(l — a/3i?HMF)^^" with a finite 
\a\ ^ 10~^ is already in complete disagreement with the 
observed dynamical fluctuations for M = 10^. We ap- 
plied the same procedure for different values of M and to 



other stationary and QSSs stemming from different ini- 
tial conditions [19j obtaining similar agreements between 
our theoretical scheme and the dynamical simulations. 

We have studied the statistical mechanics of QSSs 
emerging in the Hamiltonian dynamics of the HMF 
model in contact with a reservoir. We have shown that 
weak interparticle correlations and the CLT implies [l5| 
that the statistical mechanics in F-space is obtained by 
restricting the BG approach to a sub-manifold defined by 
a nonequilibrium value of the magnetization m — (m) j2| ■ 
During the QSS, the HMF does not thcrmalizc with the 
TB. The temperature to be used in the Boltzmann fac- 
tor is fixed by the fundamental thermodynamic relation 
applied in this nonequilibrium situation and corresponds 
to twice the specific kinetic energy of the system. Our 
theoretical approach, based on the idea of incomplete 
equilibrium given the quasi-stationary values of the 
order parameter and the temperature, allows one to cal- 
culate the other thermodynamic quantities such as the 
energy of the system and its fluctuations (i.e., the spe- 
cific heat). We expect the present approach to be signif- 
icant for nonequilibrium systems displaying stationarity 
or quasi-stationarity |^ il^ 20] concomitantly 
with a kinetic theory based on the Vlasov or Balescu- 
Lenard equations 0|. 
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